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Transport measurements on the cuprates suggest the presence of a quantum criti- 
cal point hiding underneath the superconducting dome near optimal hole doping. 
We provide numerical evidence in support of this scenario via a dynamical clus- 
ter quantum Monte Carlo study of the extended two-dimensional Hubbard model. 
Single particle quantities, such as the spectral function, the quasiparticle weight 
and the entropy, display a crossover between two distinct ground states: a Fermi 
liquid at low filling and a non-Fermi liquid with a pseudogap at high filling. Both 
states are found to cross over to a marginal Fermi-liquid state at higher tempera- 
tures. For finite next-nearest-neighbor hopping t' we find a classical critical point at 
temperature T c . This classical critical point is found to be associated with a phase 
separation transition between a compressible Mott gas and an incompressible Mott 
liquid corresponding to the Fermi liquid and the pseudogap state, respectively. 
Since the critical temperature T c extrapolates to zero as t' vanishes, we conclude 
that a quantum critical point connects the Fermi-liquid to the pseudogap region, 
and that the marginal-Fermi-liquid behavior in its vicinity is the analogous of the 
supercritical region in the liquid-gas transition. 
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1. Introduction 

(a) Relevance of quantum criticality in the cuprates 

The u nusually high supe r condu cting transition temperature of the hole doped 



cuprates (|Bednorz fc Mullen . 119861 ) remains an unsolved puzzle, despite more than 
two decades of intense theoretical and experimental research. Pairing, which has 
a d— wave symm etry and short cohere nce length, but too high of a T c to be ac- 



counted by BCS (jBardeen et all 119571 ). is not the only unconventional property of 
these materials. Their phase diagram, shown in Fig.[TJ is a landscape of exotic states 
of mat ter. Undope d cuprates are Mott insulators with antiferromagnetic long-range 



order ( Neell . 1949h . Antiferromagnetism collapses upon small doping and it is re- 



placed by a pseudogap state characterized by a suppression of spectral weight along 
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the antinodal direction. Further doping turns the system into a conventional Fermi- 
liquid metal. Between the Fermi-liquid and the pseudogap region lies a strange metal 
phase with linear-T resistivity. The superconducting dome emerges in the cross-over 
between the pseudogap and the Fermi-liquid regions at lower temperatures. 

Strong electro nic corre l ations are the cause of the rich phase diagram of cuprate 
superconductors ( Phillips . 2010() . The same strong correlations render traditional 
theoretical approaches, such as perturbation theory and Fermi-liquid theory, in- 
applicable. Some recent conceptual progress has been achieved by associating the 
optimal T c with a quantum critical point (QCP), lying underneath the su percon- 
ducting dome and co nnecting the pseudogap and the Fermi-liquid regions ( BrounL 
20081: ISachdevL l201Ch . Unlike a classical critical point, a QCP affects the behavior 
of the system in a wide range of temperatures and might explain the emergence of 
a linear-T resistivity up to room temperature. 



E xperimental evidence for a QCP comes f rom transport ([van der Marel et al 



20031: iDaou et al 
(iBernhard et all. 



2005: Plate et al 



20091 : Balakirev et al . 2009) and thermodynamic measure ments 



2001). Angle-resolved photo-emission spectros copy (ARPES) (|Shen et al 



2005) and quantum oscillation measurements (|Doiron-Levraud et al 



20071 ) show that in the pseudogap region, the Fermi surface consists of small pockets 



which have different topology than the large Fermi surface present in the Fermi liq- 
uid. It is reasonable to assume that those two states are orthogonal to one another 
and are connected through a transition. Additional evidence in support of quantum 
crit icality comes fro m measurements of the Kerr signal in YBCO by Jing Xia et 



al. (|Xia et all [2008) . They find that at the pseudogap crossover temperature, T* a 



non-zero Kerr signal develops sharply and persists even inside the superconducting 
dome . This is consi s tent w ith earlier neutron scattering measurements by Fauquc 
et aZ.I Fauaue et al . 2006f ). which show the development of magnetic order in the 



pseudogap phase. 

In this manuscript we review numerical evidence of quantum criticality in the 
Hubbard model, the de-facto model for the cuprates, that appeared in earlier pub- 
lications. In those cited works, the Hubbard model is solved using the dynamical 
cluster approximation (DCA) in conjunction with several quantum Monte Carlo 
(QMC) cluster solvers. In all calculations relevant for the phase diagram we neglect 
the superconducting transition. The interplay be tween the QCP an d superconduc- 
tivity will be discussed in a future publication. ( Yang et al . 2010l) In this review 
we focus on the thermodynamic quantities, such as the entropy and the chemical 
potential, and also on single-particle quantities, such as the spectral weight and 
the quasiparticle weight. The thermodynamic properties give unbiased evidence of 
quantum criticality, whereas single-particle properties may be used to gain more 
detailed insight on the ground state. Both set of quantities rely on the evaluation 
of the self-energy which can be calculated using quantum cluster methods. 

At a critical interaction-dependent filling, we find that the entropy exhibits 
a maximum, the quasiparticle weight displays a crossover from Fermi liquid to 
pseudogap behavior, and the spectral function shows a wide saddle point region 
crossing the chemical potential. This is consistent with the presence of a QCP, since 
the lack of an energy scale results in an enhanced entropy at low temperatures. 
We also find that by tuning an appropriate control parameter, the next-nearest- 
neighbor hoping, t' , the QCP becomes a classical critical point associated with a 
phase separation transition. We present our findings in two sections. In section [3J 
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Figure 1. The phase diagram of the 
cuprates. As a function of temperature and 
doping, the cuprates display antiferromag- 
netic order at low doping, a non-Fermi liq- 
uid pseudogap region at intermediate dop- 
ing and a metallic region at higher doping. 
Around optimal doping, superconductivity 
develops, and above the superconducting 
dome, a strange metal with non-Fermi liq- 
uid properties appears. T* separates the 
pseudogap from the marginal Fermi-liq- 
uid phase. Tx is the crossover tempera- 
ture between the Fermi and the marginal 
Fermi-liquid regions. A quantum critical 
point hides underneath the superconduct- 
ing dome near optimal hole doping. 



we discuss the single-particle spectra and the thermodynamics properties of the 
t' = Hubbard model. In section [31 we discuss the phase separation in the t' > 
Hubbard model. 

(b) Hubbard Model 

Short after the discovery of high-T c superconductors, Anderson ([Anderson 



19871 ) suggested that the Hubbard model captures t he basic proper t ies of the high 



temperature superconductors and Zhang and Rice ([Zhang fe Ricd . I1988T ) demon- 
strated that only a single band is needed. The single-band Hubbard model is rep- 
resented by the Hamiltonian: 



H = —t 



^ E [4 



H.c. 



(1.1) 



where c\ a (cj ff ) is the creation (annihilation) operator of an electron at site i and spin 
is the corresponding number operator, t is the hopping parameter between 



nearest-neighbor sites, and U the on-site Coulomb repulsion. Despite its apparent 
simplicity, the Hubbard model is not oriously difficult to solve. No analyti c al so- 
lutions exist except in on e dimension ( Lieb fc Wul . 11968k iFrahm fe KorepinL [l990; 
Kawakami fc Yane . 1990h . However, tremendous theoretical and computational ef- 
forts have resulted in approximation schemes that provide access to the physics of 
this model in higher dimensions. In this manuscript we also discuss results for the 
generalized Hubbard model which includes hopping between next-nearest neighbor 
with amplitude t': 



H 



[4 



H.c. 



E [ c l c <° 



H.c. 



(1.2) 



«M» 



Important progress in our understanding of strongly correlated models has been 
achieved by the development of finite size methods, including exact diagonalization 
and QMC. The latter works well in the simulation of bosonic systems where creation 
and annihilation operators commute. However, due to the minus sign problem asso- 
ciated with the anticommutation relations of fermionic operators, QMC is limited 
to small lattice sizes and consequently give questionable predictions for correlated 
electronic systems in the thermodynamic limit. 
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Another successful approach is the dynamical mean-field approximation (DMFA) 
which treats the local dynamical correlat ions explicitly and non-local (inter-site) 
correlations in a mean-field appr oximation ( Georges et al. . 1996; Me tzner fe Vollhardtl . 
1989; Muller-HartmannLll989QlO). This technique become s exact in the limit of in- 



finite diln^n^ioiiT^Geo^^^^^otliar , 1992t Jarrelj 1992h . However, when applied 
to finite dimensions, the DMFA fails to describe the renormalization effects due 
to momentum-dependent modes and the transitions to phases with non-local order 
parameters. Thus, DMFA misses physical phenomena that are abundant in strongly 
correlated systems, such as the development of spin or charge density wave phases, 
localization in the presence of disorder, spin-liquid physics, unconventional super- 
conductivity, etc. 

The limitations of th e DMFA are addres sed by cluster mean-field theories. Those 
fall into two cat egories dMaier et al , 20051) : the cluster dynamical mean field the- 
ory ( CDMFT) dKot iar et al. . 200lh . which is formulated in real space, and the 
DCA (jHettler et all . 19981) which is formulated in momentum space. In both cases 
the system is viewed as a cluster embedded in an effective medium. The formal 
difference between DCA and CDMFT is that in real space, the DCA cluster sat- 
isfies periodic boundary conditions whereas the CDMFT cluster is open. The two 
methods shou ld give the same re s ults f or large enough clusters. In this work we 
present DCA (jHettler et all . Il998l l2000h results. 

DCA treats short-ranged correlations explicitly, while longer ranged ones are 
approximated by the mean field. By increasing the cluster size, the length-scale 
of the explicitly treated correlations can be gradually increased while the calcula- 
tion remains in the thermodynamic limit. In momentum space, the DCA can easily 
be conceptualized as the approximation in which the self-energy calculated by the 
coarse grained g reen function. Quant um Monte Carlo based solvers such as Hirsch 
Fye ( HFQMC) (jHirsch fc Fvel . Il986l) . continuous-time (CTQ MC) dRubtsov et ~al. 



2005) and determinantal quantum Monte Carlo (DQMC) (jBlankenbecler et al. 



19811 ) are used to solve the cluster problem. QMC methods are often formulated in 



imaginary time and an analytic continuation to real time is necessary to evaluate 
physical quantitie s. Fortunately, powerful techniques such as the maxim um entropy 
method (MEM) (jGubernatis et all Il99lt Ijarrell fc Gubernatisl Il996l) are able to 
successfully select the most likely solution. 

Even though quantum cluster schemes have provided a tremendous breakthrough 
in our understanding of the Hubbard model, they are also subject to limitations. 
Quantum Monte Carlo solvers suffer from the sign problem, which scales exponen- 
tially with inverse temperature, interaction strength and cluster size. This limits the 
application of the method to relatively small cluster sizes, higher temperatures and 
intermediate interactions. The limitation in the cluster size is particularly problem- 
atic close to a phase transition where the correlation length diverges. The coarse- 
graining also limits the momentum resolution, which for typical cluster sizes is too 
small to capture detail features of the spectra, such as van Hove singularities. For 
a Fermi liquid, this is not a limitation since the physics is dominated by the low 
frequencies in which the self-energy is momentum independent. However, intrinsi- 
cally anisotropic states, such as the pseudogap, or possibly the quantum critical 
region, can be captured only approximately. Finally, MEM uses Bayesian statistics 
to find the most likely spectra for the QMC data, subject to sum rules, such as 
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conservation of the spectral weight. Because of the statistical errors in the QMC 
data, the frequency spectrum resolved using MEM has a limited resolution. 

Despite those limitations, progress can be achieved in accessing the quantum 
critical region by algorithmic optimizations. A truly universal way in dealing with 
the severity of the sign problem is to vastly increase the statistics, using massively 
parallel QMC algorithms with highly optimized codes. 



2. Prom Fermi Liquid to Pseudogap 

A great advantage of the DCA is its ability to evaluate the self-energy as a func- 
tion of momentum k and Matsubara frequency zu;„, S(k, iw n ). From the self-energy 
various single-particle quantities, such as the spectral function, A(k, lj), the quasi- 
particle weight, Zk, and the energy can be derived. All those quantities provide 
insight about the ground state of the system. In this section we will show how 
the transition from the Fermi-liquid to the pseudogap state is reflected in such 
single-particle quantities, 
(a) Spectral Function 

The single-particle spectral function shows a clear evolution from a Fermi-liquid 
to a pseudogap state as the filling increases towards half filling. Fig. [5] displays a 

density plot of the spectral function, A(k, to) = 3G(k, to), which is extracted by 

analytically continuing the imaginary time Green function. At low filling, n < 0.85, 
the spectral function exhibits a typical Fermi-liquid form. A notable characteristic 
is th e presence of a wide sad dle point region, reminiscent of a van Hove singular- 



ity, (|Radtke fc Normanl . [l99l along the antinodal direction. Around the critical 



filling of n = 0.85 this saddle point feature crosses t he chemical potential. Thi s 



crossing results in a sharp peak in the density of states ( Vidhvadhiraia et all 2009f ) 



which displays low-energy particle- hole symmetry ( Chakrabortv et all . 2008h . We 



are currently explor ing the influen ce of the van Hove singularity on the supercon- 
ducting transition. (jPathak et all ) At higher filling, n > 0.85, the spectral weight 
collapses along the antinodal direction and a pseudogap opens. The Fermi surface 
obtained by extremizing Vrik shows a similar evolution (see lower panels in Fig. 
The Fermi-liquid region consists of a large hole pocket, which extends and touches 
the edges of the Brillouin zone (0, ±7r), (±7r, 0) at n — 0.85. In the pseudogap region 
the Fermi surface consists of four Fermi arcs centered around the nodal points, sim- 
ilar to the ones obtained from ARPES. These results clearly demonstrate that the 
DCA can capture qualitatively the evolution of the ground state from a Fermi-liquid 
to a pseudogap phase. 

(b) Quasiparticle Weight 

Whereas the spectral function gives a qualitative understanding of the ground 
state, it relies on the analytic continuation of numerical data. Since extracting 
quantitative information from analytically continued data is difficult, a more robust 
way is to rely on imaginary time quantities, such as the quasiparticle weight Z(k). 
Since the quasiparticle weight is finite across a Fermi surface, but vanishes if the 
spectrum is incoherent, it will allows to clearly distinguish between a Fermi liquid 
and a pseudogap state. The quasiparticle weight can be directly obtained from the 

1 - I , where wo = 

w J 
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Figure 2. Upper panels: Density plots of the spectral function A(k,uj) for the Fermi 
liquid (left), marginal Fermi liquid (middle) and pseudogap region (right) for filling 
n = 0.75, 0.85 and 0.95, respectively. (The dashed feature seen in the regions of steep- 
est dispersion, especially for n = 0.75, is a plotting artifact). The momentum is along the 
path G(0, 0) — > M(7T, 7r) — >• X(n, 0) — > G(0, 0). A wide saddle point region between X and 
G sits above the chemical potential in the Fermi-liquid region and crosses it around the 
critical filling (n = 0.85). In the pseudogap region this features sits below the chemical 
potential leaving a gap along the antinodal direction behind it. Note that the fact that the 
dispersion looks discontinuous along G(0,0) — > M(n,n) in the left and middle panels is 
an artifact of our interpolation algorithm. Lower panels: Fermi surface as extracted from 
|Vnk| in the Fermi liquid (left), marginal Fermi liquid (middle) and pseudogap (right) re- 
gion showing the development of the pseudogap in the antinodal direction. The Coulomb 
repulsion is U = 6t, the temperature T = 0.069t, and the cluster size N c = 16. The energy 
unit is it. 



is the lowest fermionic Matsubara frequency. At the limit T — >• and for a well- 
behaved self-energy, Zq (k) converges to the quasiparticle weight, Z (k). Fig. [3] (a) 
displays Zan — Zq(ujo — nT,k \\ (0,0) — > (0, 7r)), the Matsubara quasiparticle 
weight alo ng the antinodal momentu m direction for U = 6t and a cluster of size 
N c = 16 (|Vidhvadhiraia et all [2009). The momentum k at the Fermi surface is 
determined by maximizing |Vn(k)|. Zan exhibits two distinguishable behaviors: 
for n > n c — 0.85 the quasiparticle weight vanishes, whereas it approaches a finite 
value for n < n c . The n > n c region corresponds to the pseudogap state in which 
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Figure 3. a) The antinodal quasiparticle fraction Zan as a function of temperature for 
different values of filling, U = 6t and cluster size N c = 16 (the unit of energy is At). The 
onset of the pseudogap region is determined by the vanishing of the antinodal spectral 
weight at zero temperature. The dashed and solid lines represent fits of the low tem- 
perature (T < 0.3) data to marginal Fermi liquid (red solid curves), Fermi liquid (black 
solid curves) and crossover forms (dashed black curves), respectively. The arrows show 
the corresponding crossover temperatures Tx and T* . The value of T* presented here is 
obtained from the spin susceptibility as explained in l|Vidhvadhiraia et aZ.L [2009T ) . but is 
consistent with the one extracted from the from the fitting forms. The ratio Zm/Zan 
of the quasiparticle weight in the nodal ((n,n)) and antinodal ((0, tt)) directions (inset) 
diverges as the pseudogap develops in accordance with Fig. [5] b) The crossover tempera- 
tures Tx and T* as a function of filling as extracted from the temperature dependence of 
Zan l|Vidhvadhiraia et all 120091 1 for the same parameters. 



the spectral weight collapses along the antinodal direction, while the n < n c region 
behaves as a Fermi liquid. 

The temperature dependence of Zan (Fig- [SI (a)) not only provides information 
about the ground state but also allows the extraction of relevant energy scales. 
By comparing the numerical results with analytical expressions derived from par- 
ticular phenomenological forms of the self-energy, we obtain Tx and T*. At low 
filling, n < n c , the high T dependence of Zan is best fit by a marginal Fermi-liquid 
form, whereas for low T the data is best fit by a Fermi liquid. The crossover oc- 
curs at a temperature Tx, which is extracted by fitting with a crossover function, 
and is accompanied by a change in the sign of the curvature of Zan- At higher 
filling (n > 0.85), the high temperature Zan can also be fit by a marginal Fermi 
liquid, whereas at low temperatures, it cannot. The crossover temperature T* can 
be extracted as the lowest temperature where the marginal Fermi liquid fit lies 
within the statistical error. However a more accurate value can be obtained from 
the bulk spin susceptibi lity which exhibits a peak at T* and the two values are 



peak ai 

found to be consistent (jVidhvadhiraia et all l2009h . The crossover temperatures 



Tx and T* are shown in Fig. [3] (b) . Both of them converge to zero as the filling 
approa ches n c = 0.85, which is the same value for which the peak in the density of 
states (jVidhvadhiraia et ~al\. l2009h crosses the chemical potential, 
(c) Thermodynamics 

A different perspective at the transition from a Fermi liquid to the pseudogap 
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Figure 4. a) The filling dependence of the entropy divided by temperature, S/T, for various 
temperatures at U = 6t and N c = 16. With decreasing temperature a peak develops around 
the critical filling of n c = 0.85. b) The temperature dependence of the chemical potential 
(i for different fillings. At the critical filling, n c , fj, becomes temperature independent at 
low temperatures. 

state comes from the evaluation of the entropy. We obtain the entropy by integrating 
the energy using the formula: 

rP 

S{(3,n) = S(0,n)+pE{(3,n)- E{0 ,n)dp , (2.1) 

Jo 

where ft is inverse temperature and S(0,n) is the infinite temperature entropy. 
Eauation l2.1l is appropriate for QMC calculations, because the integration reduces 
the statistical error. The challenge is to have good en ough statistics to contr ol the 
error of the surface term, (3E{fi, n). In Mikelsons et al. (jMikelsons et a/.l . l2009t ) large 
statistics was possible simply by using large computational resources. The entropy 
divided by the temperature, shown in Fig. 0] (a) , exhibits a maximum at exactly the 
same critical filling that was identified before from the spectral function and the 
quasiparticle weight. In Fig. 01(b), we show the chemical potential, function 
of temperature. We note that at the critical filling dfi/dT = 0, since the entropy 
and the chemical potential are related by the Maxwell relation: 

(-) =-(^) ■ (2.2) 
\dn) TU \dT) Un 1 } 

Also the temperature dependence of the chemical potential can be used as a 
practical criterion to identify the location of the critical filling, because evaluat- 
ing the chemical potential is much less computationally intensive than evaluating 
the entropy. Using this criterion we investigate the important question of the de- 
pendence of n c on the Coulomb repulsion U. As it is shown in Fig. [5J we find 
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Figure 5. a) The chemical potential as a function of temperature for fillings of n = 0.85 
and 0.90 and for a variety of interaction strengths U for N c = 12. b) The critical filling, 
defined by the filling in which du/dT = versus U. The critical filling decreases with U 
monotonically and is projected to reach the atomic limit value of n c = 2/3 at U c = 30t. 



that increasing U reduces the critical filling and thus enlarges the pseudogap re- 
gion in the phase diagram . Our results follow the trend proposed in earlier argu- 
ments ( Chakrabortv et al . 2008) according to which the critical filling decreases in 
order to reach the atomic limit value of n c = 2/3. 

In this section we have shown that several single-particle quantities are con- 
sistent with the presence of a QCP. The qualitative form of the single-particle 
spectrum shown in Fig. [5] is fundamentally different in the Fermi-liquid and the 
pseudogap regions, which points to orthogonal ground states. The temperature 
dependence of the quasiparticle weight reveals the presence of two crossover tem- 
peratures T* and Tx, which converge to zero at n c as shown in Fig. [3] (b). If the 
crossover temperatures Tx and T* constitute energy scales that suppress degrees 
of freedom, their vanishing at n c means that there are no relevant energy scales to 
quench the entropy and therefore it collapses at a slower rate, which is consistent 
with the peak of the entropy observed at n c . The natural next step to investigate 
quantum criticality is to access the QCP. However the fermion sign problem severely 
limits the applicability of quantum Monte Carlo techniques close to a QCP. It is 
possible however, as we will discuss in the next section, that by tuning an appro- 
priate control parameter, the critical point may be lifted to finite temperature and 
thus studied with QMC. 



3. Phase Separation and Quantum Criticality 

Experiments suggest that cuprate superconductors are susc eptible to charge in - 



homogeneities, such as stripes or checkerboard modulations (jHinkov et all 12004) . 
These inhomogeneous charge patterns have stimulated intense theoretical and ex- 
perimental research. Here we will consider the possibility that those charge insta- 
bilities are evidence that the cuprates are close to a phase separation transition, 
and this proximity may be related to the nature of the QCP. 

Our findings suggest that the Hubbard model displays a phase diagram similar 
to the one for the gas-liquid transition with Mott liquid (ML) and Mott Gas (MG) 
regions. Fig. Ela) shows a possible phase diagram for the Hubbard model as a 
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Figure 6. a) The schematic phase diagram in the presence of charge separation. This phase 
diagram describes the transition between two states labeled Mott Liquid (ML) and Mott 
Gas (MG) as a function of temperature, T, chemical potential, fi, and filling, n. The red 
surface represents the coexistence region, which terminates in a critical point (CP). As 
we go around the critical point the state changes smoothly from ML to MG. Along the 
first order transition line and for a fixed T and u, the filling has two values, b) Filling 
as a function of chemical potential for several temperatures in the vicinity of the charge 
separation critical point. The number next to each curve represents the temperature. The 
coexisting phases are an incompressible Mott liquid at n « 1 and a compressible Mott 
gas at n ~ 0.93. The critical temperature is T c — O.lt. The blue dashed line represents 
the surface of metastability which is not accessible within the DCA. The green dotted line 
represents the isothermal of the metastable state inside the phase coexistence region (gray 
zone). At the critical point the isothermals for T > T c cross. The inset shows the scaling 
curve (n — n c ){T — T C )~P vs (/i — fi c )(T — T c )~ l3S in arbitrary units for fi c = 3i, n c = 0.96, 
Tc — O.lt. The scaling exponents, ft = 0.10 ±0.05 and ft 8 ~ 1, are roughly consistent with 
the Ising universality class. 



function of T, and n. The red-colored surface is a schematic of the region where 
the Mott liquid and Mott gas states, characterized by different densities, coexist 
for T < T c . The critical point is located at temperature T c , filling n c , and chemical 
potential /i c . One can go from one state to the other either smoothly, by avoiding 
the phase separation region, or through a first-order transition by crossing it. Right 
on the phase separation region, the density has two values for a given value of /i 
and T. 

Macridin et al. (jMacridin et all . [2006) provided compelling evidence of phase 
separation in the case of the generalized Hubbard model (Eq. (|1.2p ) with positive 
next-near-neighbor hopping t' — 0.3 1 and U = 8t. Using the DCA in a N c = 8 
cluster with HFQMC as the cluster solver, they showed that below a critical tem- 
perature T c ~ O.lt a first order transition occurs, which is identified by a hysteresis 
in the n versus /i curve for T < T c . As shown in Fig. H>tb) with more precise data 
obtained using DQMC as the cluster solver, the hysteresis is between two states 
of different filling, the Mott liquid at half filling and the Mott gas at a filling of 
about 0.93 for T = 0.07U. The Mott liquid is incompressible and insulating. Its 
compressibility, which is the slope of the filling vs \i curve in the high filling side 
of the hysteresis curve, is small and decreases with temperature. Also the density 
of states of the ML phase, shown in Fig. [TJa), exhibits a gap as expected for an 
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Figure 7. The density of states of the a) Mott liquid and b) Mott gas states at T = 0.077t 
(dotted line) and T = 0.057 (solid line). The Mott liquid is an incompressible insulator 
with a pseudogap while the Mott gas is weakly compressible with a Fermi liquid peak in 



insulator. On the other hand, the Mott gas is compressible and metallic; the density 
of states is finite at the chemical potential (fj, = u> = 0), as displayed in Fig. E£b). 

The analogy to the well-known phase diagram of a liquid-gas mixture, such as 
water and steam, is useful to understand this phase transition. At low temperatures, 
there is a region in the pressure-volume phase diagram in which water and steam 
coexist for a range of pressures. As the temperature is increased, the region of coex- 
istence contracts and finally terminates at a critical point where the compressibility 
diverges. In the pressure-temperature phase diagram, this region of coexistence be- 
comes a line of first order transitions which terminates at a second order point where 
the water and gas become indistinguishable and the compressibility diverges. Since 
the line terminates, it is possible for the system to evolve adiabatically from steam 
to water without crossing a phase transition line; therefore, the steam and water 
must have the same symmetry. 

In the Mott liquid and Mott gas system the chemical potential /i replaces the 
pressure and the density n replaces the volume of the water-gas mixture. Because 
the order parameter separating the ML from the MG, the density n, does not 
have a continuous symmetry, order may occur at finite temperatures, and the ML- 
MG transition will most likely be in the Ising or lattice gas universality class. 
Within this context, one may then understand the hysteresis of Fig. [B{b). The solid 
lines are isotherms which show how the system evolves with increasing density. 
At the temperature T — T Cl the compressibility diverges at the critical filling. 
As the temperature is lowered further, there is a region where the ML and MG 
coexist. Inside this region the isothermals contain unphysical regions of negative 
compressibility (dashed green line in Fig. [6fb)) along with metastable regions of 
positive compressibility. The metastable branch of the isothermal in the vicinity 
of the ML is a "supercooled" ML, whereas the one in the vicinity of the MG is a 
"superheated" MG. The translational invariance of DCA along with the stabilizing 
effect of the mean-field host enable access to those metastable states. However 
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Figure 8. (a) The chemical potential-temperature phase diagram of the ML and MG 
mixture for t' > 0. The ML and MG coexist on a line of first order transitions with 
positive slope. Since ML and MG have the same symmetry this line can terminate in a 
second order critical point. The blue dashed lines define the boundaries of the supercritical 
region where the ML and MG cannot be distinguished. Outside this region either the ML 
or MG character dominates, (b) The chemical potential-temperature phase diagram for 
i! — i 0. The first-order line is absent but supercritical region remains as a quantum critical 
region. In the Hubbard model the lines T* and Tx (Fig. EJb)) define the boundaries of 
this region. 



the real physical system will phase separate and the two phases will co-exist in 
equilibrium (dotted blue line in Fig. Efb)) . 

We can sketch the phase diagram in the T — fi plane using the analogy to the 
water-steam mixture. The most generally applicable rule governing the shape of 
phase diagrams was established by Gibbs. For a system of c conserved components 
and / phases, the Gibbs constraint is give by the relation $ = c — / + 2 where $ 
is the number of independent variables needed to specify the state of every phase. 
In this case, as in the water-steam system, the number of components c = 1, since 
the particle number is conserved. At a location in the phase diagram where only 
one phase exists, $ = 1 — 1 + 2 = 2, so there are two independent variables, and 
the phase diagram is a surface on the /i, T and n three-dimensional space. There 
will be places in the phase diagram where two phases exist simultaneously, then 
<I> = 1 — 2 + 2 = 1, implying that two phases co-exist only along lines in the 
phase diagram. At the lines in the T — fi plane where two phases coexist, n is also 
determined for each phase, but its value can be different. That is a line of first order 
transitions. 

Additional information about the lines of first order transitions is obtained from 
the Clapeyron's equation. The Gibbs free energy G = E — TS — fiN , and dG = 
—SdT — Ndji, must be the same for the coexisting phases on a line. If we label the 
two phases 1 and 2, then 

(Si - S 2 )dT = -{Ni - N 2 )dfi . (3.1) 

If we identify the latent heat L = (Si — S2)T, then dfi/dT = —L/(TAn) represents 
the slope of the first order transition line. Since the latent heat L of going from 
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Figure 9. (a) Filling, n, vs. chemical potential, for T = 0.077t, N c — 16, U = 6t and 

various is shown in solid lines and the compressibility — in dashed lines. A critical 

filling, identified by the peak in the compressibility appears at higher temperatures and 
fillings as t' is increased. The inset shows the t! dependence of the critical filling, n c . (b) 
Schematic phase diagram of the Hubbard model in the ^t, t' and T space (neglecting su- 
perconductivity). The classical critical point turns asymptotically into a quantum critical 
point as t' — > 0. 



ML to MG is positive, but dn is negative, the slope d/i/dT of the line of first order 
transitions is positive. 

Above the critical point terminating the ML-MG transition, the system displays 
supercritical behavior in a region where the gas and the liquid cannot be distin- 
guished thermodynamically (c.f. Fig. [5]). It is possible for the system to evolve adi- 
abatically through a counterclockwise path from deep in the MG region, through 
the supercritical region, into the ML region. At the lower edge of the supercritical 
region, the system loses the Fermi-liquid character of the MG, and at the upper 
edge, it begins to acquire the pseudogap character of the ML. 

Let us discuss now how this phase separation, which occurs at finite temperature, 
is related to quantum criticality. The key parameter is the next-nearest-neighbor 
hoping, t' . For t' = there is no evidence for phase separatio n at finite T, but such 



a phase separation occurs for positive t' . Khatami et al. ( Khatami et al. . 2010h 
performed a systematic analysis of the phase diagram of the extended Hubbard 
model as a function of t' . As shown in Fig. [9] (a) the compressibility, \c — dn/dfi. 
exhibits a peak for all positive t' at a critical filling that depends on t' . The width of 
the peak measures the distance from the critical temperature: the sharper the peak 
the closer to T c the employed temperature is. We see that the critical temperature 
increases with t' and it starts from T c = at t' = 0. These results point to the 
phase diagram of Fig. [5] (b). At a positive t' a charge separation occurs at tem- 
peratures T < T c (t') and at a critical filling n c (t') between an incompressible and 
insulating Mott liquid and a compressible metallic Mott gas. Right at T c , there is a 
terminating second-order critical point. By decreasing t' this second-order critical 
point is pushed down to lower temperatures. Presumably the line of second-order 
critical points terminates at the QCP. 

Such a scenario constitutes a new path to quantum criticality as it is closely 
tied to charge fluctuations rather than spin fluctuations. However, numerous sim- 
ulations suggest that a finite positive t' enhances antiferromagnetic correlations, 
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and since phase separation is only present for t' > it suggests that it is driven by 
strong spin correlations. In addition, previous simulations incorporating Holstein 
phonons to the H ubbard model found that phonons also enhance the phase sepa- 
ration instability ( Macridin fc Jarrell . 20091 ). As t' jt — > (and the electron-phonon 
coupling vanishes) , the phase separation critical point approaches zero temperature 
becoming a QCP. Here, the first-order behavior is absent from the phase diagram 
(Fig. [3f b)) leaving only the adiabatic path from the ML to the MG, which passes 
through the supercritical region, which is now the quantum critical (QC) region. 
The crossover scale Tx and the pseudogap scale T* are now understood as the 
boundaries of the QC region. As we cross the line of T* from the QC region into 
the ML region, the characteristics of the ML become apparent, including the pseu- 
dogap in the DOS and the insulating behavior. As we cross the line of Tx from the 
QC region into the MG, the characteristics of the MG become apparent, including 
Fermi liquid formation. 

Those calculations certainly do not elucidate the nature of the Mott liquid and 
Mott gas states in real materials. The long-ranged nature of the Coulomb inter- 
action prevents true charge separated states, but the phase separation we observe 
may also correspond to other charge instabilities, such as stripes or checkerboard 
patterns. To distinguish between different charge instabilities, systematic calcula- 
tions in much larger clusters are necessary which are not practical at the moment. 
However, whatever the type of order, those calculations provide convincing evidence 
for the existence of a first-order transition at low temperatures. Such a transition 
is similar to the liquid-gas or the ferromagnetic transition and its phase diagram 
would look like Fig. E£a): a first order line of coexistence which terminates at a 
critical point at a critical temperature T c and a critical filling n c . 



4. Conclusions 

The presence of a QCP at finite filling in the cuprate phase diagram is a topic of 
active theoretical and experimental research. Quantum cluster methods are able 
to shed some light in this phase diagram. By studying single-particle quantities 
for t' = 0, such as the spectral function and the entropy, it can be shown that a 
Fermi-liquid region at low filling and the pseudogap region at higher filling have 
different spectral signatures, and are connected through an intermediate "marginal 
Fermi- liquid" region of maximal entropy. Due to limitations of quantum Monte 
Carlo, the ground state and quantum criticality are not accessible. We also neglect 
the superconducting phase transition. The connection with quantum criticality is 
established by switching on t'. For positive t' a classical critical point emerges at 
finite temperature T C) which increases with t'. We note that t' is not the only control 
parameter that may be able to tune the critical point to finite temperatures, but 
other parameters, such as phonon coupling, may have the same effect. The phase 
diagram around the critical point is similar to that of the gas-liquid transition, where 
the incompressible Mott liquid and the compressible Mott gas are the coexisting 
phases. The strange metal region in this context may be viewed as the supercritical 
region lying in the vicinity of the critical point. Within the scenario we presented, 
the pseudogap region is not characterized by an order parameter, rather it must 
have the same symmetry as the Fermi-liquid and the marginal Fermi-liquid, since 
these regions are connected by an adiabatic path in the T — fi phase diagram. 
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Further investigation is necessary to fully characterize the pseudogap region, and 
also to investigate the connection of those results with other scenarios of quantum 
criticality. 
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